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1. Introduction 

Our understanding of the phase structure of QCD at high baryon density and low temperature 
remains severely hampered by the sign problem. In the absence of first-principles methods which 
have been proven to circumvent this problem, we can study a related theory, QCD with colour 
group SU(2) (QC2D), which does not suffer from the sign problem. This may firstly allow us to 
confront model studies with lattice results, thereby constraining these models in their application to 
real QCD, and secondly reveal generic features of the phase structure of strongly interacting gauge 
theories, including the nature of deconfinement at high density. 

Here we present an update of our ongoing investigation of the phase structure of QC2D as a 
function of temperature and chemical potential [jl], ^. 

2. Simulation details 

We study two-colour QCD with a conventional Wilson action for the gauge fields and two 
flavours of unimproved Wilson fermion. The fermion action is augmented by a gauge- and iso- 
singlet diquark source term which serves the dual purpose of lifting the low-lying eigenvalues of 
the Dirac operator and allowing a controlled study of diquark condensation. Further details about 
the action and the Hybrid Monte Carlo algorithm used can be found in [^. We have performed 
simulations at j8 = 1.9, k: = 0.168, corresponding to a lattice spacing a = 0.178fm, determined from 
the string tension, and a pion mass aniji = 0.645 or niji « VlOMeV. The ratio of the ground state 
pseudoscalar to vector masses is mj^/mp = 0.80 [Q]. Our lattice volumes and the corresponding 
values for temperatures T, chemical potentials n and diquark sources j are given in table [l[ All 
results shown are extrapolated to 7 = using a linear Ansatz except where otherwise stated. 







T (MeV) 


Ha 


ja 


16 


24 


47 


0.3-0.9 


0.04 


12 


24 


47 


0.25-1.1 


0.02, 0.03, 0.04 


12 


16 


70 


0.3-0.9 


0.04 


16 


12 


94 


0.2-0.9 


0.02, 0.04 


16 


8 


141 


0.1-0.9 


0.02, 0.04 



Table 1: Lattice volumes and associated temperatures T, chemical potentials and diquark sources j. 



3. Order parameters and phase structure 

The left panel of fig. |l| shows the diquark condensate (qq) = {y^"'CY5T2V^ — V^CY5Z2¥^'^) 
a function of chemical potential, for the A^^ = 24, 12 and 8 lattices. In the case of a weakly-coupled 
BCS condensate at the Fermi surface, the diquark condensate, which is the number density of 
Cooper pairs, should be proportional to the area of the Fermi surface, ie {qq) ~ jJ.^. 

For the lowest temperature, T = 41 MeV (A^^ = 24), we see that (qq) /}X^ has a plateau in the 
region 0.35 < jtifl' < 0.6. The increase for jia > 0.6 may be evidence of a transition to a new state of 
matter at high density, although the impact of lattice artefacts cannot be excluded. The lower Umit 
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Figure 1: Left: The diquark condensate {qq) I pi^ extrapolated to j = Q for = 24, 12,8 (r = 47,94, 141 
MeV). Right: The renormalised Polyakov loop as a function of chemical potential, for all temperatures. The 
shaded symbols are for ja — 0.04; the open symbols are ja = 0.02. The filled black circles are the results 
for the 16^ x 24 lattice. The dashed line indicates the inflection point of L at jU = 0. The inset shows the 
unrenormalised Polyakov loop. 



of the plateau roughly coincides with the onset chemical potential jx^ ^ nij^/l ^ 0.33a^', below 
which both the quark number density and diquark condensate are expected to be zero. We find no 
substantial volume dependence at any jU. Our results at 7 = 70 MeV (not shown here) are almost 
identical to those at T = 47 MeV. At T = 94 MeV (A^^- = 12), (qq) is significantly suppressed, and 
drops dramatically for jxa > 0.7. At T = 141 MeV (A^^- = 8) the diquark condensate is zero at all 
/I, confirming that the system is in the normal phase at this temperature. 

In the right panel of fig. [I] we show the order parameter for deconfinement, the Polyakov loop 
L, for our four temperatures. It has been renormalised by requiring L{Ta = = 0) = 1, see 
for details. We see that for each temperature T, L increases rapidly from zero above a chemical 
potential pid{T) which we may identify with the chemical potential for deconfinement. In the 
absence of a more rigorous criterion, we have taken the point where L crosses the value it takes at 
Tii{n = 0), Ld = 0.6 to define }JLd{T)- The results are shown in fig. ^ with error bars denoting 
the range Ld =0.5-0.7. To more accurately locate the deconfinement line, we will need to perform 
a temperature scan for fixed jU -values. This is underway. 

The estimates of critical chemical potentials for deconfinement and superfluidity can be trans- 
lated into a tentative phase diagram, shown in fig. In summary, from the order parameters we find 
signatures of three different regions (or phases): a normal (hadronic) phase with {qq) = 0, (L) ^ 0; 
a BCS (quarkyonic) region with {qq) ~ at low T and intermediate to large /i; and a deconfined, 
normal phase with {qq) = 0, (L) 7^ at large T and/or pL. After extrapolating our results to 7 = we 
see no evidence of a BEC region described by ^PT, with {qq) ^ \/\— jifjix'^ iQ], in contrast with 
earlier work with staggered lattice fermions |^]. This may be related to the large value of mji/mp 
in this study. Simulations with lighter quarks may yield further insight into this. 

In the right panel of fig. |2| we show the static quark potential computed from the Wilson loop 
at Nt: = 24, for /la = 0.3,0.5,0.7,0.9. We find that as we enter the superfluid region, the string 
tension is slightly reduced, but that this is reversed as n is increased further, leading to a strongly 
enhanced string tension at jxa = 0.9, which according to our analysis of the Polyakov loops should 
be in the deconfined region. This agrees with the pattern that was already observed in [^. We also 
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Figure 2: Left: A tentative phase diagram, including the location of the deconfinement transition in the 
(/i, r) plane, determined from the renormalised Polyakov loop, and the transition to the diquark condensed 
{qq) ^ phase. Right: The static quark potential computed from the Wilson loop, for the 12^ x 24 lattice 
and different values of ji. 
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Figure 3: The quark number density divided by the density for a noninteracting gas of lattice quarks (left) 
and continuum quarks (right). 



find no significance 7-dependence in our results. At present we do not have a good understanding 
of why the static quark potential should become antiscreened at large /i. Computing the static 
quark potential using Polyakov loop correlators rather than Wilson loops may yield further insight 
into this issue. 



4. Equation of state 

We now turn to the bulk thermodynamics of the system, and in particular the quark number 
rig and the energy density e. Fig. ^ shows the quark number density rig for A^^ = 24, 12 and 8, 
extrapolated to zero diquark source, and normalised by the noninteracting value for lattice fermions 
on the left and for continuum fermions on the right. The difference between the two gives an 
indication of the lattice artefacts. We see that the density rises from zero at jj. ^ Ho = 0.32a^' , and 
for the two lower temperatures is roughly constant and approximately equal to the noninteracting 
fermion density in the region 0.4 < /xa < 0.7. The peak at /la ~ 0.4 in the A^^- = 24 data in the upper 
panel is an artefact of the normalisation with hsb for a finite lattice volume: the raw numbers for 
the 12-' X 24 and 16-' x 24 lattices are identical within errors, but usb differs by about 50% around 
jxa = 0.4. 
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Figure 4: The quark and gluon contributions to the energy density (left) and total energy density (right), 
divided by /x^, for ja = 0.04 (open symbols) and j — (filled symbols). 



The density for = i does not show any plateau as a function of /x; instead, riq/nsB shows 
a roughly linear increase in the region 0.4 < /la < 0.7. This is suggestive of the system being 
in a different phase at this temperature. We also note that riq/nsB for Nt = 12 rises above the 
corresponding A^^ = 24 data for jxa > 0.7, where, according to the results of Sec. ||, the hotter 
system is entering the deconfined, normal phase. 

These results lend further support to our previous conjecture that in the intermediate-density, 
low-temperature region the system is in a "quarkyonic" phase: a confined phase (all excitations are 
colourless) that can be described by quark degrees of freedom. 

The renormalised energy density can be derived by going to an anisotropic lattice formulation 
with bare anisotropics jg = y^^t/l^s, 7q = ^t/ Ks and physical anisotropy ^ = Qs/a^. In the isotropic 
Umit 7^ = 7g = i§ = 1 the energy density is then given by £ = £g + £^ with 
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K 
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We have determined the Karsch coefficients dcj/ dt, with c,- = 7^ , 7^ , j8 , K" by performing simulations 
with 7^, 7g / 1. Our estimates for these coefficients are 



0.90-14, 



+40 

0.13-5, 



dl 



+0.24 

0.59-1.37, 



dK_ 



+69 

-0.052-15. 



(4.3) 



Our results for the energy density are shown in fig. ^. We see that the quark contribution is 
negative for all values of /i and T , but this is balanced by the positive gluon contribution, giving 
a positive or zero total energy. The energy density is very sensitive to the values of the Karsch 
coefficients [|l[]; for example, if dyq/dt, is changed from the suprisingly low value of 0.13 to a 
more 'natural' value of 0.8, we find that £^ > for jxa > 0.6. 



5. Heavy quarkonium 

We have investigated the heavy quarkonium spectrum by computing non-relativistic QC2D 



correlators on our A^^ = 24, 16 and 12 lattices. We use an 



lattice NRQCD lagrangian 
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Figure 5: Left: Temperature dependence of the '5o state energy vs. jU for Ma = 5.0 with j = 0.04. Right: 
The ratio Y.xG{fi;T; ji) / Y,,^G{x, t;0) for ^Pq correlators on 12^ x 24 with Ma = 5.0. Due to the noisiness of 
the P-wave data, only a limited T range is shown. 



to compute the heavy quark Green function; see for further details. We find that the S-wave 
correlators can be fitted with an exponential decay oc ^-^^"^ even once /i 7^ 0; moreover the fits 
are quite stable over large ranges of T, indicating that 5-wave bound states persist throughout the 
region 47 MeV < T < 90 MeV. 

Fig. ^ shows the T- and /i -dependences of the state energy AE. We see that as is 
varied, initially the state energy decreases from that at jj. =0, but once /i reaches the region 

0.5) < }xa < /i2(— 0.85), the state energy stays roughly constant. For }x > }X2, the 
state energy starts increasing again. In contrast to the observables studied in Sees |3| and ^ we find 
no clear, systematic dependence on the diquark source term for jxa < 0.5. For jxa > 0.5 on the 
other hand, AE{ja = 0.02) < AE{ja = 0.04). This suggests that the energy, extrapolated to j = 0, 
may continue to decrease up to jxa « 0.7 before increasing. 

As the temperature increases from 47 MeV (A^^ = 24) to 70 MeV (A^^ = 16) we find that 
the point where the energy of the state starts increasing goes from /ic? ?a 0.7 to 0.55. This is 
consistent with the estimate of the deconfinement transition in Sec. ^ For A^^ = 12 we do not yet 
have any data in the -region which might confirm this. It is interesting to note that AE increases 
with increasing T, in accordance with what has been observed in hot QCD with = 

In contrast to the S'-waves, it is difficult to find stable exponential fits to the P-wave correlators 
with the current Monte-Carlo data before statistical noise sets in, except for the case jxa < 0.25. In 
the right panel of fig. || we instead show the ratios of the ^Po correlators at different values fJ. y^O 
to the correlator at = 0. Note that any effect we observe is entirely due to the dense medium. 

The S-wave correlator ratios show an increase with T which corresponds to the negative ^^o 
energy difference A£'(/x) — A£'(/x — 0) that was previously observed. In the quarkyonic region, the 
P-wave ratios behave similarly to the 5-wave, but in the deconfined region (n > /I2), the P-wave 
ratios are non-monotonic, initially decreasing with T before turning to rise above unity for r/a 4. 
On the other hand, the P-wave correlator ratios on the 12^ x 16 and 16^ x 12 lattice show monotonic 
behavior similar to those of the S-waves, suggesting a subtle interplay of density and temperature 
effects on the P-wave states. 
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6. Summary and outlook 

From lattice simulations of dense QC2D at a range of temperatures, we have identified three 
distinct regions of the phase diagram: a hadron gas at low /i and T, a quarkyonic region at inter- 
mediate n and low to intermediate T, and a deconfined quark-gluon plasma at high T and/or jj.. 
Taking the limit of zero diquark source has served to make our identification of the quarkyonic 
region more robust. Investigations into the exact nature and location of the deconfinement and the 
superfluid to normal transitions are underway, as are simulations at smaller lattice spacings and 
with smaller quark masses. 
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